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Abstract 

A series of experiments measuring high-energy cosmic rays have recently reported strong indica- 
tions for the existence of an excess of high-energy electrons and positrons. If interpreted in terms 
of the decay of dark matter particles, the PAMELA measurements of the positron fraction and 
the Fermi LAT measurements of the total electron-plus-positron flux restrict the possible decaying 
dark matter scenarios to a few cases. Analyzing different decay channels in a model-independent 
manner, and adopting a conventional diffusive reacceleration model for the background fluxes of 
electrons and positrons, we identify some promising scenarios of dark matter decay and calcu- 
late the predictions for the diffuse extragalactic gamma-ray Hux, including the contributions from 
inverse Compton scattering with the interstellar radiation field. 
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I. INTRODUCTION 



Different experiments measuring high-energy cosmic rays have over the last months re- 
ported a wealth of new results pointing to the existence of an exotic source of electrons and 
positrons. The PAMELA collaboration reported evidence for a sharp rise of the positron frac- 
tion at energies 7 — 100 GeV possibly extending toward even higher energies, compared 
to the expectations from spallation of primary cosmic rays on the interstellar medium 2|. 
This result confirmed previous hints about the existence of a positron excess from HEAT 3|] , 
CAPRICE and AMS-01 Jsj. Almost at the same time, the balloon-borne experiments 
ATIC 6| and PPB-BETS 7| reported the discovery of a peak in the total electron-plus- 
positron flux at energies 600 — 700 GeV, while the H.E.S.S. collaboration js] reported a 
substantial steepening in the high-energy electron-plus-positron spectrum above 600 GeV 
compared to lower energies. 

These results raised a lot of interest in the astrophysics and particle physics communities, 
leading to many proposals trying to explain this excess. One of the most popular astrophys- 
ical interpretations of the positron excess is in terms of the electron-positron pairsproduced 
by the interactions of high-energy photons in the strong magnetic field of pulsars jo, 1^, 11 1. 



However, this interpretation requires a rather large fraction of the spin-down power being 
injected in the form of electron-positron pairs or a rather large rate of gamma-ray pulsar 
formation. Alternatively, the positrons could be originating from the decay of charged pions, 
which are in turn produced by the hadronic interactions of high-energy protons accelerated 
by nearby sources [l^ . 

An arguably more exciting explanation of the cosmic-ray positron excess is the possibility 
that the positrons are produced in the annihilation or the decay of dark matter particles. 
Should this interpretation be confirmed by future experiments, then the positron excess 
would constitute the first non-gravitational evidence for the existence of dark matter in our 
Galaxy. The interpretation of the PAMELA excess in terms of dark matter is subject to 
constraints from the flux measurements of other cosmic-ray species. A very important con- 
straint arises from the measurements of the antiproton flux by PAMELA 13|, BESS95 |14| . 
BESS95/97 [15], CAPRICE94 [16], CAPRICE98 [ItI and IMAX 



which are consistent 



with the expectations from conventional propagation models, thus excluding the possibility 
of a large antiproton flux from dark matter annihilation or decay 19|, 1201 ]. 
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The steep rise in the positron fraction observed by PAMELA can be explained by dark 
matter annihilations in the center of the Milky Way, provided the dark matter particle has a 
mass larger than ~ 150 GeV and annihilates preferentially into leptons of the first or second 
generation ,21|. This interpretation of the positron excess, however, typically requires the 
ad hoc introduction of large boost factors. Furthermore, it has been argued that if dark 
matter annihilations are the origin of the PAMELA anomaly, then the predicted gamma- 
ray emission from the center of the Galaxy is in conflict with the H.E.S.S. observations for 
typical cuspy halo profiles |22|. On the other hand, if the positron excess is due to the decay 
of dark matter particles, the dark matter particles must have a mass larger than ~ 300 
GeV, a lifetime around 10^^ s, and must decay preferentially into hard leptons of the first 



or second generation 



231]. In this case, no boost factors are required and the gamma and 



24]. Some recent works on 
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28]. 



radio measurements are consistent with present measurements 
the indirect detection of decaying dark matter can be found in 

More recently, the Fermi LAT collaboration has published measurements of the electron- 
plus-positron flux from 20 GeV to 1 TeV of unprecedented accuracy 29[ , revealing an energy 
spectrum that roughly follows a power law oc E^'^'^ without any prominent spectral fea- 
tures. Simultaneously, the H.E.S.S. collaboration reported a measurement of the cosmic-ray 
electron-plus-positron spectrum at energies larger than 340 GeV, confirming the Fermi result 
of a power-law spectrum with spectral index of 3.0 ± 0.1(stat.) ± 0.3(syst.), which further- 
more steepens at about 1 TeV |30] . The measured energy spectrum is harder than expected 
from conventional diffusive models, although it can be accommodated by an appropriate 
change of the injection spectrum of primary electrons. However, when taken together with 
the steep rise in the positron fraction as seen by PAMELA up to energies of 100 GeV, the 
Fermi LAT data suggest the existence of additional Galactic sources of high-energy electrons 
and positrons with energies up to a few TeV. Furthermore, it should be borne in mind that 
the determination of the correct Galactic cosmic-ray scenario is still an open problem, and 
while an electron injection spectrum harder than the conventional could reproduce the Fermi 
data, it fails to account for the AMS-01 and HEAT data below 20 GeV and the H.E.S.S. 
data above 1 TeV |3l[ ]. 



In this paper we analyze the constraints that the results of the PAMELA and Fermi 
collaborations impose on the scenario of decaying dark matter, assuming a GALPROP 
conventional model as our Galactic cosmic-ray scenario. To this end, we pursue a model- 
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independent approach, calculating the prediction for the positron fraction and the total 
electron-plus-positron flux for various decay channels of both a fermionic and a bosonic dark 
matter particle. We will identify the most promising scenarios in the light of the PAMELA 
and Fermi data, and we will calculate for those the predictions for the antiproton flux and the 



diffuse extragalactic gamma-ray flux. Some related works have recently appeared |32l. |33|. 

The paper is organized as follows: in Section 2 we will review the production and propa- 
gation in the Galaxy of high-energy electrons/positrons, antiprotons and gamma rays from 
dark matter decay, including a contribution to the total gamma-ray flux from inverse Comp- 
ton radiation. In Section 3 we will show the predictions for the positron fraction and the 
total electron-plus-positron flux for several decaying dark matter scenarios. For the promis- 
ing scenarios, we will also show the predictions for the antiproton and the gamma-ray fluxes. 
Finally, in Section 4 we will present our conclusions. 

II. COSMIC RAYS 

In this section, we briefly review the propagation model for cosmic rays that we need 
for the calculation of the electron, positron and antiproton fluxes measurable at Earth. 
Furthermore, we discuss our calculation of the flux of gamma rays, which come from inverse 
Compton scattering (ICS) with the interstellar radiation field (ISRF) as well as directly from 
the decay process itself. 

If dark matter decays at a sufficiently large rate, the decay products (electrons, positrons, 
antiprotons and gamma rays) could be observable as an anomalous contribution to the high- 
energy cosmic-ray fluxes. The production rate of particles per unit energy and unit volume 
at a position r with respect to the center of the Milky Way is given by 

Q(Ef)= — (1) 

where dN/dE is the energy spectrum of particles produced in the decay and p(r) is the 
density profile of dark matter particles in the Milky Way halo. For definiteness we will 
adopt the spherically symmetric Navarro-Frenk- White halo density profile [3^: 

(r/rc)[l + (r/rc)J^ 

with Po — 0.26GeV/cm^ and ~ 20 kpc, although our results are almost independent of 
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choice of the density profile^. 



A. Electron/positron propagation 

After being produced in the Milky Way halo, the electrons and positrons propagate 
through the Galaxy and its diffusive halo in a rather complicated way before reaching the 
Earth. The propagation is commonl y d escribed by a stationary two-zone diffusion model 
with cylindrical boundary conditions 35|] . Under this approximation, the number density of 
electrons and positrons per unit energy, fg±{E,r,t), satisfies the following transport equa- 
tion: 

= ^ = V ■ [K{E,r)Vf,±] + ^[b{E,f)f,^] + Q,±(E,f) . (3) 

The first term on the right-hand side of the transport equation is the diffusion term, which 
accounts for the propagation through the tangled Galactic magnetic field. The diffusion co- 
efficient K{E, r) is assumed to be constant throughout the diffusion zone and is parametrized 
by K[E) = Kq 13 TZ^, where (3 = v/c and IZ is the rigidity of the particle, which is defined 
as the momentum in GeV per unit charge, TZ = p{GeV)/Z. The second term accounts for 
energy losses due to ICS on starlight or the cosmic microwave background (CMB), syn- 
chrotron radiation and ionization. We parameterize the energy loss rate as b{E) = -^i^, 
with Eq = 1 GeV and te = 10^^ s. Lastly, Qe±{E,r) is the source term of electrons and 
positrons, defined in Eq. ([T]). 

The boundary conditions for the transport equation, Eq.Q, require the solution 
/e±(-E', r, t) to vanish at the boundary of the diffusion zone, which is approximated by a 
cylinder with half-height L = 1 — 15 kpc and radius R = 20 kpc. Under these assumptions, 
the propagation of electrons and positrons can be described by just three parameters, the 
normalization Kq and the spectral index 6 of the diffusion coefficient, which are related to the 
properties of the interstellar medium, and the height of the diffusion zone, L. In our numer- 
ical analysis we will adopt for these parameters the values of the MED propagation model 



defined in 



361 ]. which provide the best fit to the Boron-to-Carbon (B/C) ratio: 6 = 0.70, 



^ Due to the effective energy loss of electrons, the high-energy component of the spectrum mostly originates 
from sources within the Galactic neighborhood of a few kpc from the Solar System, where the different 
halo profiles are very similar. We have checked that choosing different halo profiles has a negligible effect 



on our results (see also [25[). 
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Kq = 0.0112 kpc^/Myr and L = 4kpc. Our conclusions, however, are rather insensitive to 
the choice of propagation parameters, as the different sets of propagation parameters yield 
rather similar results for cosmic rays from local sources. This is due to the fact that at 
higher energies above several 10 GeV energy losses dominate the effects of diffusion, render- 
ing the exact propagation model parameters less relevant. We illustrate the dependence of 
the results on the adopted model parameters for a particular example in Section 3. 

The solution of the transport equation at the heliospheric boundary, r = Tq, z = 0, can 
be formally expressed by the convolution 

/e± (E) = dE' G,. {E, E') ^^M^ . (4) 

Mdm tdm Jo dE' 

The Green's function Ge±{E, E') encodes all the information about astrophysics (such as the 

details of the halo profile and the propagation of electrons/positrons in the Galaxy), while 

the remaining part is model-dependent and is determined by the nature of the dark matter 

particles. Analytical and numerical expressions for the Green's function for the propagation 



of electrons/positrons can be found in 25 1. 



Finally, the interstellar flux of primary electrons/positrons from dark matter decay is 
given by: 

^e^E) = ^fe^{E) . (5) 

In order to compare our results with the PAMELA results of the positron fraction as well 
as the Fermi results on the total flux of electrons plus positrons it is necessary to know the 
background fluxes of high-energy electrons and positrons. The background flux of positrons 
is constituted by secondary positrons produced in the collision of primary protons and other 
nuclei with the interstellar medium. On the other hand, the background flux of electrons 
is constituted by a primary component, presumably produced in supernova remnants, as 
well as a secondary component, produced by spallation of cosmic rays on the interstellar 
medium and which is much smaller than the primary component. Whereas the spectrum 
of secondary electrons and positrons is calculable in a given propagation model, the energy 
spectrum and the normalization of the primary electron flux is unknown and has to be 
determined by direct measurements. 

In this paper we will adopt for the background fluxes of electrons and positrons the ones 
corresponding to the "model 0" presented by the Fermi collaboration in [3l|, which fits well 
the low-energy data points of the total electron-plus-positron flux and the positron fraction. 



The interstellar background fluxes can be parametrized as: 

^'(g)- ( ,f^;;;,3 )GeV--m-'s--s.--. (6) 

- ( ifo'.ooll^^ + 2" GeV- sr- . (7) 

where e = E/1 GeV. In the energy regime between 2 GeV and 1 TeV these approximations 
are better than 5%. 

At energies smaller than ~ 10 GeV the electron/positron fluxes at the top of the atmo- 
sphere can differ considerably from the interstellar fluxes, due to solar modulation effects. 



Under the force fleld approximation 
to the interstellar fluxes by the following simple relation 38l |: 



37|, the fluxes at the top of the atmosphere are related 



<^I^HEtoa) = ^<^'ME,s), (8) 

where £"13 = -Etoa + 4>f, with £"13 and -Ktoa being the electron/positron energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively, and being 
the solar modulation parameter, which varies between 500 MV and 1.3 GV over the eleven- 
year solar cycle. In order__to compare our predictions with the AMS-01 and HEAT data we 
will take (pp = 550 MV 

Then, if there exists an exotic source of electrons and positrons from dark matter de- 
cay, the total flux of electrons plus positrons and the positron fraction at the top of the 
atmosphere read, respectively, 

= $°^(£) + $f+^(£) + k <l>l'l\E) + <I>'^+HE) , (9) 
PF(E) = ^ , f ^ ^ 10) 

where we have left the normalization of the primary electron flux as a free parameter, k, 
to be determined in order to provide a qualitatively good flt to the PAMELA and Fermi 
measurements. 



B. Antiproton propagation 

Antiproton propagation in the Galaxy can be described in a similar manner as that 
of electrons and positrons. However, since antiprotons are much heavier than electrons 
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Model 


6 


Ko (kpcVMyr) L (kpc) (km/s) 


MIN 


0.85 


0.0016 


1 13.5 


MED 


0.70 


0.0112 


4 12 


MAX 


0.46 


0.0765 


15 5 



TABLE I: Astrophysical parameters compatible with the B/C ratio that yield the minimal (MIN), 
median (MED) and maximal (MAX) flux of antiprotons. 



and positrons, energy losses are negligible. However, antiproton propagation is affected by 
convection, which accounts for the drift of antiprotons away from the disk induced by the 



Milky Way's Galactic wind. Following j36| we will assume that it has axial direction and 
that it is constant inside the diffusion region: Vc(r) = K sign(2;) k. Then, the transport 
equation for antiprotons reads: 

df- 

= ^ = V • [K{T, f)Vfp - V,{f)fp] + Qp{T, f) , (11) 

where T is the antiproton kinetic energy. 

As for the case of electrons and positrons, the solution of the transport equation at the 
heliospheric boundary, r = r©, z = 0, can be formally expressed by the convolution 

UiT) = jr^—— / dT' G,{T, T') , (12) 

where T^ax = ^dm ~ and mp is the proton mass. Analytical and numerical expressions 



for the Green's function Gp{T,T') can be found in 25|]. Finally, the interstellar flux of 



primary antiprotons from dark matter decay is given by 

$r(^) = i^fp(T) ' (13) 

where v is the velocity of the antiprotons. The prediction of the antiproton flux from dark 
matter decay is very sensitive to the choice of propagation parameters. Therefore, we will 
show the results for three different propagation models that are consistent with the observed 
B/C ratio and that give the maximal (MAX), median (MED) and minimal (MIN) antiproton 
flux [36]. The relevant parameters are summarized in Tab. [H 

The antiproton flux at Earth is also affected at low energies by solar modulation effects. 
Again, under the force fleld approximation js?!, the antiproton flux at the top of the Earth's 
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atmosphere is related to the interstellar antiproton flux 38|] by the simple relation: 

where Tis = Ttoa + 4>f, with Tjs and Ttoa being the antiproton kinetic energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively. 



C. Gamma rays from inverse Compton scattering 

As discussed above, electrons and positrons from dark matter decay lose their energy 
mainly via interaction with the Galactic magnetic fleld and the ISRF. In the flrst case 
(assuming injection energies of the order of 1 TeV) synchrotron radiation in the radio band 
with frequencies 0(0.1 — 100 GHz) is produced and potentially observable (see e.g. Ref. 
In the second case, the ICS of electrons and positrons on the ISRF (which includes the 
cosmic microwave background, thermal dust radiation and starlight) produces gamma rays 
with energies between 100 MeV and 1 TeV. Recently, ICS in connection with the PAMELA 



39|) 



excess was discussed in Refs. 



4Q|, 



review of ICS can be found in Ref. 



IL and we will follow their treatment. A pedagogical 
42]. 



The rate of inverse Compton scattering of an electron with energy Ee, where an ISRF 
photon with an energy between e and e + de is upscattered to energies between E.y and 
+ dE^, is given by 



dN{Ee,f) 3 dT 



/iSRF(e, r) 



1 (gF) 



2 



2glng + 1 + g - 2g2 + ^-^^(1 - q) 

21 + ql 



(15) 



de dE-y dt 4 7^ e 

where ut = 0.67 barn denotes the Compton scattering cross section in the Thomson limit, 
7 = Ee/rrie, q = E^rrie/ {'ie'~^{me'~^ — E^)), and /isRF(e, '^) denotes the differential number 
density of ISRF photons with energy e, at spatial position r. 

From Eq. (|T5l) one can derive the energy loss bics{Ee,r) of electrons due to ICS (which 
represents the dominant contribution to the energy loss rate b{Ee, r) in the transport equa- 
tion, Eq. ([3])), and the corresponding power V{E^, Ee,r) that is emitted in gamma rays 
with energies between E^ and E^ + dE^ via 

hcs{E^, r) = r de P dE.iE, - e)'^4^^ (16) 

and 



^ ~e de dEy dt 



V(E^,E.,r>^fME^-^/^^. (17) 



Neglecting diffusion and synchrotron losses^, the energy distribution of electrons and 
positrons from dark matter decay it is given by 



DM 



/e±(^e,r) = j-^r-— / dE Q,±{E,r). (18) 

For the differential flux of ICS photons of energy E^ from a region AQ of the sky it then 
follows: 

"^-^ 1 f ,af r- y(E.l (19) 



dE^ Ait E^tbm J An Ao.s. M-dm Jm, bics{Ee,r) 

where Y{Ef.) = dE dNf,± / dE describes the number of particles in the spectrum of 

electrons and positrons above a certain energy E^. In the second integral, the coordinate s 
runs over the line of sight (l.o.s.), which points in the direction of Afi. The prefactor 2 takes 
into account that the same amount of gamma rays comes from the dark matter electrons 



and positrons. In this work, we use the ISRF data as derived in Ref. [4^], and we fully take 
into account the spatial dependence of the energy loss bics{Ee,r) in Eq. (fTSi) . Furthermore, 
we calculate the gamma rays from ICS with extragalactic origin. In this case, effects of 
redshifting must also be taken into account. Details of this calculation can be found in 



Ref. [41 1 . 



In addition to the gamma-ray signal from dark matter decay there exists a background 
contribution, presumably originating from active galactic nuclei (AGN), which is perfectly 
isotropic, and which has an energy spectrum which is assumed to follow a simple power law; 
the normalization and index will be treated as free parameters to be determined by requiring 
a good fit of the total flux to the data. 

We will compare our predicted flux to two sets of data for the extragalactic diffuse gamma- 
ray background obtained from the EGRET data, using two different models for the Galactic 
background, averaging over the whole sky, excluding the region of the Galactic plane with 
latitudes \ b\ < 10°. The first analysis of the extragalactic diffuse gamma-ray flux by Sreeku- 
mar et al. 



4J] revealed a power law 



dE 



J bkg 



1.37 X 10"^ (g^) s)"' (20) 



In Ref. '32'| it was discussed that this approximation gives results for the ICS fluxes that are correct at 
the 0(2) level, which is sufficient for our analysis. 
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in the energy range 50 MeV — 10 GeV. On the other hand, the extraction of the extragalactic 



background by Strong, Moskalenko and Reimer |45|, using an optimized model to better 



simulate the Galactic diffuse emission, revealed a steeper power law. 



J Tl / 7-. \ -0.32 

aJ 



dE 



6.8 X 10"^ ( ^^^^ ^^'^ ^^^^ 



'J bkg 

between 50 MeV — 2 GeV and an intriguing break of the spectrum at energies 2 GeV — 
10 GeV. Future measurements by the Fermi LAT, as well as a better understanding of the 
Galactic diffuse emission, will provide a better determination of the extragalactic diffuse 
gamma-ray flux in the near future. 

III. PREDICTIONS FROM DECAYING DARK MATTER 

We will analyze in this section the predictions for the positron fraction and the total 
electron-plus-positron flux including a possible contribution from dark matter decay in or- 
der to account for the anomalies observed by PAMELA and Fermi. To keep the analysis 
as model-independent as possible, we will analyze several scenarios of decaying dark mat- 
ter, computing the predictions for the positron fraction and the total electron-plus-positron 
flux for either a fermionic or a bosonic particle, which decays into various channels with a 
branching ratio of 100%. We calculated for each of these channels the energy spectrum of 
electrons and positrons using the event generator PYTHIA 6.4 4^. Thus, from the particle 
physics point of view the only free parameters are the dark matter mass and lifetime. From 
the astrophysics point of view there are a number of uncertainties, such as the choice of 
propagation parameters and the choice of the background fluxes of electrons and positrons. 
As mentioned in the previous section, we will adopt the MED propagation model defined in 



361], which provides the best fit to the Boron-to-Carbon (B/C) ratio, although the results 
are not very sensitive to the particular choice of the propagation model. On the other hand, 
for the background fluxes of electrons and positrons we will adopt the spectra corresponding 
to the "model 0" proposed by the Fermi collaboration. However, we will allow for a possible 
shift in the normalization of the background flux of electrons, which is dominated by pri- 
maries, due to our ignorance of the amount of electrons injected in the interstellar medium. 
In our analysis we will sample several dark matter masses and treat the dark matter lifetime 
and the normalization of the background flux of electrons as free parameters which will be 
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determined to provide a qualitatively good fit to the PAMELA and Fermi measurements. 
Note that below energies of 10 GeV the data is best fitted for normalizations /c ~ 1. In our 
plots, we always used normalization factors k > 0.8. 

Let us now discuss the cases of fermionic and scalar dark matter particles separately. 



A. Fermionic dark matter decay 

In the case where the dark matter particle is a fermion V'dm, we consider the following 
decay channels^: 



V 



(22) 



where the three-body decay into charged leptons and a neutrino is assumed to be mediated 
by the exchange of a scalar particle, motivated by the interesting scenario of a hidden gaugino 



as dark matter particle 



28|. 
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FIG. 1: Positron fraction {left panel) and total electron-plus-positron flux (right panel) for the 
decay channel Vdm Z^i^ with Mdm = 100 TeV (solid) and 5TeV (dotted). The dashed line 
shows the background fluxes as discussed in the text. Solar modulation is taken into account using 
the force field approximation with (pp = 550 MV. 

The predicted positron fraction in the case where the dark matter particles decay via 
^DM Z'^i^ is shown in the left panel of Fig. [H compared to the PAMELA, HEAT, 



^ We do not include quarks or Higgs bosons in the list, since they yield similar signatures to gauge boson 
fragmentation. Furthermore, we only consider decay channels with two or three final-state particles. 
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FIG. 2: Same as Fig. [H but for the decay channels V'dm — *■ W^£^. Upper panels: V'dm W^e^ 
with My)m = 2000 GeV (sohd) and 300 GeV (dotted). Middle panels: V'dm ^ W^n'f with Mdm = 
3000 GeV (sohd) and 600 GeV (dotted). Lower panels: Vdm ^ W^t^ with Mdm = 8000 GeV 
(sohd) and 1000 GeV (dotted). 

CAPRICE and AMS-01 data, for the exemplary dark matter masses Mdm = 5 and 100 TeV. 
In the right panel we show the corresponding total electron-plus-positron flux compared to 
the results from Fermi, H.E.S.S., PPB-BETS, BETS, ATIC, HEAT, CAPRICE and AMS- 
01. The dark matter lifetimes and the normalization factors k of the primary electron flux 
have been chosen in each case to provide a reasonable flt to the PAMELA and Fermi data 
points. In this decay channel, the dominant source of electrons and positrons is the frag- 
mentation of the boson (with a rather small branching ratio of decays into a pair 
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FIG. 3: Extragalactic diffuse gamma-ray ffux for V'DM W^fi^ with Mdm = 3000 GeV and 
TDM = 2.1 X 10^^ s. The gamma-ray flux is averaged over the whole sky, excluding the Galactic 
plane, \b\ < 10°. We included gamma rays produced directly in the final-state radiation of the 
muons and the fragmentation of (green line), gamma rays from ICS of dark matter electrons 
and positrons on the ISRF (solid blue line; the dotted blue lines show, from left to right, the fluxes 
that come from scattering on the CMB, on the thermal radiation of dust and on starlight) and 
gamma rays from ICS outside of our Galaxy (red). The black solid line shows the overall flux. The 
dark red and dark blue lines show the total flux (dash-dotted) adding an isotropic extragalactic 
background (dashed) with a power-law spectrum. Normalization and power index are chosen to fit 



one of the two data sets shown 



of charged leptons), which produces relatively soft particles. As a result, even though this 
decay mode can produce a visible excess in the positron fraction, the energy spectrum is in 
general too flat to explain the steep rise observed by PAMELA. An exception occurs if the 
dark matter mass is very large, M^m ^ 50 TeV. In this case, the electrons and positrons 
from dark matter decay are boosted to high enough energies to produce the steep rise in the 
positron fraction. However, these large dark matter masses seem to be in conflict with the 
H.E.S.S. observations, which require a falloff in the total electron-plus-positron spectrum at 
~ 1 TeV. 

On the other hand, we show in Fig. [2] the predictions for the cosmic-ray electron and 
positron fluxes when a fermionic dark matter particle decays as ipDU W'^i'^ for different 
dark matter masses. The electrons and positrons created in the fragmentation of the 
gauge bosons produce a rather fiat contribution to the positron fraction. However, the hard 
electrons and positrons resulting from the decay of the /i^ and leptons or directly from 
the dark matter decay into produce a rise in the total energy spectrum and in the positron 
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FIG. 4: Antiproton flux (left panel) and the corresponding antiproton-to-proton ratio {right panel) 
for V'DM W^fj,^ witli Mdm = 3000 GeV and tdm = 2.1 x 10^^ s. For the antiproton flux we adopt 



the background from Ref. 



from Ref. 



47l |. while antiproton-to-proton ratio is plotted using the background 



48| . and the yellow band indicates the uncertainties from the propagation model. The 



solid black line corresponds to the MED model of Tab. [H 
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FIG. 5: Same as Fig. [H but for the flavor-democratic decay ipBM — * W^£^ with equal branching 
ratios into the three charged lepton flavors, for Mdm = 2000 GeV (solid) and 300 GeV (dotted). 

fraction. The decay mode tpBM —>■ W^e^, which can produce a steep rise in the positron 
fraction and is thus consistent with the PAMELA observations, produces also a steep rise 
and a sharp falloff in the total electron-plus-positron flux, which is not observed by Fermi. 
Thus, the possibility that the dark matter particles decay preferentially in this decay mode, 
which is well compatible the PAMELA observations, is now strongly disfavored by the Fermi 
results on the total electron-plus-positron flux. 

The decay mode ipBU W^fi^, however, can nicely accommodate the PAMELA and 
Fermi observations when the dark matter mass is Mdm — 3 TeV and the lifetime is tdm — 
2.1 X 10^^ s. In this decay mode, the fragmentation of the gauge bosons also produces 
fluxes of primary antiprotons and gamma rays, which are severely constrained by present 
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FIG. 6: Same as Fig. [H but for the decay channels V'dm i^i^iy. Upper panels: iI^bm — ^ e e'^u 
with Mdm = 2000 GeV (sohd) and 400 GeV (dotted). Middle panels: V'dm fJ.^fJ.^iy with Mdm = 
3500 GeV (sohd) and 1000 GeV (dotted). Lower panels: Vdm ^ r"r+i/ with Mdm = 5000 GeV 
(sohd) and 2500 GeV (dotted). 

experiments. The predictions for the gamma-ray and antiproton fluxes for this particular 
decay mode are shown in Figs. [3] and HI the former figure shows the gamma-ray fluxes 
from final-state radiation and fragmentation (green), and from Galactic (blue) and 
extragalactic (red) ICS of dark matter electrons and positrons. We also show the total flux 
cornpared to the extraction of the extragalactic diffuse gamma-ray flux by Sreekumar et 



al. 



44l | and by Strong, Moskalenko and Reimer 45], averaging over the whole sky excluding 



the region of the Galactic plane with latitudes \b\ < 10° and assuming a power law for 
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FIG. 7: Same as Fig. [U but for the democratic decay V'dm ~^ l^F^v with equal branching ratios 
into the three charged lepton flavors, with Mdm = 600 GeV (dotted) and 2500 GeV (solid). 
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FIG. 8: Same as Fig. [3l but for Vdm /U fi^v {left panel, with Mdm = 3500 GeV) and for the 
democratic decay V'dm l^l'^f {right panel, with Mdm = 2500 GeV). 

the genuinely extragalactic component. On the other hand, the latter figure shows the 
prediction for the antiproton-to-proton ratio with an uncertainty band corresponding to the 
MAX, MED and MIN models in Tab. [H While the absolute flux is compatible with existing 
measurements, it is apparent from the figure that the antiproton-to-proton ratio is in some 
tension with the results at the highest energies explored by PAMELA. The fragmentation 
of the gauge bosons also produces a sizable contribution to the total gamma-ray flux at 
high energies which could be visible by the Fermi LAT as a bump over the background, which 
is assumed to follow a simple power law, especially if it has a large index as in the extraction 
of the diffuse extragalactic background from the EGRET data by Strong, Moskalenko and 
Reimer. Lastly, the decay mode V'dm — ^ W^t^ predicts, for a wide range of dark matter 
masses, a positron fraction and an electron-plus-positron flux that are too fiat to explain 
the anomalies observed by PAMELA and Fermi. 

In some decaying dark matter scenarios, the dark matter particles decay into charged 
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leptons of different flavors and not exclusively in just one channel. As an illustration of the 
predictions of this class of scenarios, we show in Fig. [5] the positron fraction and the total 
electron-plus-positron flux for a dark matter particle that decays democratically into the 
three flavors, for Mdm = 2000 GeV (sohd) and 300 GeV (dotted). Although these scenarios 
could explain the PAMELA excess, the predicted spectral shape of the total flux is not 
consistent with the Fermi data: either the energy spectrum falls off at too low energies 
or it presents a sharp peak at high energies, due to the large branching ratio into hard 
electrons and positrons. Scenarios with smaller branching ratio into electron flavor and larger 
branching ratio into muon flavor could, however, explain both anomalies simultaneously. 

The dark matter particles could also decay into three fermions, namely into a lepton- 
antilepton pair and a neutrino. In this case many possibilities could arise depending on the 
specific particle physics scenario. We will just concentrate on the case where the lepton and 
the antilepton carry the same flavor and the decay is mediated by a heavy scalar^. The 
results for the positron fraction and the total electron-plus-positron flux are shown in Fig. [6l 

The spectrum produced in the decay into electron-positron pairs is flatter in this case than 
in the two-body decay ipBU VV^e^, although it still predicts a rather prominent bump in 
the electron spectrum at high energies, which is not observed by Fermi. More promising is 
the decay channel ipBM — ^ /i"//"*"//, which can reproduce quite nicely the Fermi electron-plus- 
positron spectrum and the steep rise in the positron fraction observed by PAMELA when 
the dark matter mass is Mdm — 3500 GeV and the lifetime is tdm — 1.1 x 10^^ s. Lastly, 
decays into tau flavor can qualitatively reproduce the steep rise in the positron fraction for 
dark matter masses above ~ 2.5 TeV, although as apparent from Fig. [6l lower right panel, 
the resulting electron-plus-positron spectrum has an energy dependence much steeper than 
E~^'^ at high energies, in tension with the Fermi measurements. In this case an additional 
source of high-energy positrons, coming e.g. from pulsars, must be invoked in order to 
reproduce the Fermi energy spectrum. 

As for the two-body decays ipBU W^t^, the dark matter particle could also decay 
into charged fermions with different flavor. We illustrate such a situation showing in Fig. [7] 
the predictions for the positron fraction and the total electron-plus-positron flux when the 
dark matter particles decay democratically into the three flavors, for dark matter masses 

Our results are not very sensitive to the mass splitting between dark matter particle and virtual scalar. 
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M-DM = 600 GeV (dotted) and 2500 GeV (solid). In particular, this is the case in scenarios 
where dark matter neutralinos decay into light hidden gauginos via kinetic mixing, or vice 
versa 28|]. It is interesting that these scenarios can simultaneously explain the PAMELA 



and Fermi anomalies when the dark matter mass is Mqm — 2500 GeV. For the two cases of 
three-body decay into muon flavors and democratic decay, we show the predictions for the 
extragalactic diffuse gamma-ray fluxes in Fig. [HI In both cases, they are consistent with the 
present data and show a deviation from the putative power-law behavior of the astrophysical 
background, which could be observed by the Fermi LAT, depending on the precise spectrum 
of the genuinely extragalactic contribution to the flux. 

We summarize our results for the promising fermionic dark matter scenarios, together 
with the corresponding dark matter masses and lifetimes, in Tab.HTl The impact of choosing 
other sets of propagation parameters is illustrated in Fig. [9] for the decay mode ipBM — *■ 
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FIG. 9: Illustration of the dependence on the choice of transport parameters. Same as Fig. [6l 
middle panels, but only for a dark matter mass of 3500 GeV. The solid, dashed and dotted lines 
correspond to the MED, MAX and MIN model parameters, respectively. The results for the MED 
and MAX model are very similar because the height of the diffusion zone becomes irrelevant above 
a few kpc for high-energy electrons from local sources. 
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B. Scalar dark matter decay 



For a scalar dark matter particle, we will discuss the following decay channels^: 



'DM - 




'DM - 


^ w+w 


'DM - 





(23) 

We show in Figs. [TU] and [TT] the positron fraction and the total electron-plus-positron flux 
for a scalar dark matter particle that decays exclusively into weak gauge bosons 0dm ~^ Z^Z^ 
and 0DM — * W^W~ for dark matter masses Mdm = 2TeV and lOTeV. As generically 
expected from decays into weak gauge bosons, the electrons and positrons produced are 
relatively soft, resulting in a positron fraction which is too fiat to explain the steep rise in 
the spectrum observed by PAMELA (for an exception with large dark matter masses see 
Fig. [I]). 

Decays into harder electrons and positrons can arise in scenarios where the scalar dark 
matter particle decays into a lepton-antilepton pair. We show in Fig. [12] the predictions 
for the positron fraction and the total electron-plus-positron flux when the scalar dark 
matter particle decays into fermions of the same generation, for dark matter masses between 
^DM = 300 GeV and 5 TeV. The decay 0dm e+e" can explain the steep rise in the 
positron fraction observed by PAMELA. However, it is apparent from Fig. [12] that the 
dark matter decay into this channel cannot be the origin of the Fermi excess in the total 
electron-plus-positron flux. The situation is similar if one considers democratic decay into 
all three flavors as shown in Fig. [13] Decays into softer electrons and positrons, as in the 
case when the dark matter particles decay exclusively via 0dm A*'*'/^") are more promising. 
In particular, a scalar dark matter particle with a mass Mdm — 2500 GeV and a lifetime 
tdm — 1-8 X 10^^ s, which decays exclusively into jj,'^ ji'^ pairs, can reproduce both the steep 
rise in the spectrum observed by PAMELA and the total electron-plus-positron spectrum 
measured by Fermi. The same holds true for decay into tau flavors, with Mdm — 5000 GeV 
and Tdm — 0.9 x 10^^ s. For these two decay channels, we also show the predictions for the 
gamma-ray fluxes in Fig. [TU which are again compatible with the present data and present 



^ Again, we do not include quarks or Higgs bosons in the list. Three-body decay modes like (/)dm — * 7 
are expected to give results similar to the fermionic dark matter case. 

20 



a spectral shape which could be visibly different from a power law, depending on the index 
of the genuinely extragalactic contribution. 

A summary of our results can by found in Tab. [ITl Note that one of the largest uncertain- 
ties that enter in the determination of the dark matter lifetime comes from the determination 



of the local dark matter density (see Ref. 49j for a recent analysis) since this quantity is 
inversely proportional to the corresponding flux of cosmic rays. 
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FIG. 10: Same as Fig. [U but for the decay channel (j)BM —>■ Z^Z^ with Mdm = lOTeV (solid) and 
2TeV (dotted). 
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FIG. 11: Same as Fig. [H but for the decay channel c^dm ^ W+W with Mdm = lOTeV (solid) 
and 2TeV (dotted). 



IV. CONCLUSIONS 

In some well-motivated dark matter scenarios, the dark matter particles are unstable 
and decay with a lifetime much longer than the age of the Universe. In this paper we 
have investigated whether the anomalies in the positron fraction and the total electron-plus- 
positron flux reported by the PAMELA and the Fermi LAT collaborations, respectively, 
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FIG. 12: Same as Fig. [H but for the decay channels (pBM (■'^^ ■ Upper panels: 4)dm e"'"e 
with Mdm = 2000 GeV (solid) and 300 GeV (dotted). Middle panels: (/)dm ^ /"^/^^ with Mdm = 
2500 GeV (sohd) and 600 GeV (dotted). Lower panels: (?;>dm ^ r+r" with Mdm = 5000 GeV 
(soHd) and 2000 GeV (dotted). 



could be interpreted as a signature of the decay of dark matter particles. We have shown 
that some decaying dark matter scenarios can indeed reproduce the energy spectra of the 
positron fraction and the total flux reasonably well, while being at the same time consistent 
with present measurements of the antiproton flux and the diffuse extragalactic gamma-ray 
flux. The most promising decay channels for a fermionic or a scalar dark matter particle are 
listed in Tab. [Tll where we also show the approximate mass and lifetime which provide the 
best fit to the data. It should be borne in mind that the astrophysical uncertainties in the 
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FIG. 13: Same as Fig. [H but for the decay channel (pBM ~* > democratic decay into three 
charged lepton flavors, with Mdm = 2000 GeV (sohd) and 300 GeV (dotted). 
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FIG. 14: Same as Fig. [3l but for cpuM fJ-^fJ- {left panel, "with Mdm = 2500 GeV) and ^divI 
T+T~ {right panel, with Mdm = 5000 GeV). 

propagation of cosmic rays and in the determination of the background fluxes of electrons 
and positrons are still large. Besides, the existence of a possibly large primary component 
of electrons/positrons from astrophysical sources, such as pulsars, cannot be precluded. 
Therefore, the precise values of the dark matter parameters can vary. The present results 
can nevertheless be used as a guidance for building models with decaying dark matter as an 
explanation of the PAMELA and Fermi anomalies. 

Future measurements of the extragalactic diffuse gamma-ray flux by the Fermi LAT will 
provide important information about the decaying dark matter scenario. First, since the 
Earth is located far from the center of the Milky Way halo, an anisotropy in the diffuse 
extragalactic gamma-ray flux is expected which could be observed by the Fermi LAT 26|. 
Moreover, all scenarios in Tab. HTl predict a departure from a simple power law in the energy 
spectrum of the extragalactic diffuse background, the deviation depending on the spectrum of 



the genuinely extragalactic contribution originating presumably from AGN. The observation 
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Decay Channel 


Mdm [GeV] 


TDM [1026s] 




3500 


1.1 




2500 


1.5 


0DM fJ'^fJ'^ 


2500 


1.8 


(/'DM T+r" 


5000 


0.9 


V'DM W^IX^ 


3000 


2.1 



TABLE II: Decay channels for fermionic and scalar dark matter, V'DM and 0dm j respectively, that 
best fit the Fermi and PAMELA data for the MED propagation model and the NEW halo profile. 
As discussed above, the dependence on the halo profile is negligible, while the dependence on the 
adopted propagation parameters is illustrated in Eig. [9] for the decay V'dm — > ^i^^^u. The decay 
mode ipT>M W^fj,^ is in tension with the PAMELA results on the antiproton-to-proton ratio, as 
mentioned in the text. 

of such a deviation would provide support for the decaying dark matter scenario and may 
help to discriminate among the different possibilities in Tab. UTI 
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